Collective traffic-like movement of ants on a trail: 
dynamical phases and phase transitions 
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The traffic-like collective movement of ants on a trail can be described by a stochastic cellular 
automaton model. We have earlier investigated its unusual flow-density relation by using various 
mean field approximations and computer simulations. In this paper, we study the model following 
an alternative approach based on the analogy with the zero range process, which is one of the few 
known exactly solvable stochastic dynamical models. We show that our theory can quantitatively 
account for the unusual non-monotonic dependence of the average speed of the ants on their density 
for finite lattices with periodic boundary conditions. Moreover, we argue that the model exhibits a 
continuous phase transition at the critial density only in a limiting case. Furthermore, we investigate 
the phase diagram of the model by replacing the periodic boundary conditions by open boundary 
conditions. 

PACS numbers: 



> 



O 



-a 
c 

o 
o 



X 



I. INTRODUCTION 

It is well known that no phase transition takes place 
in one-dimensional systems in thermodynamic equilib- 
rium if the interactions are short-ranged. But, systems 
of interacting particles driven far from equilibrium can 
exhibit transitions from one dynamical phase to another 
even in one-dimensional space with only short-ranged in- 
teractions P, Qj- A class of models, which have been 
receiving lot of attention in the recent years from this 
perspective, consist of interacting particles each of which 
is self-propelled rather than driven by external field 0; 
the dynamics of these systems are most often formulated 
in terms of updating rules and, therefore, these models 
can also be regarded as stochastic cellular automata (CA) 
ill 

Several different mechanisms are known to give rise 
to the phase transitions in systems of driven interacting 
particles. The boundary-induced phase transitions 6] 
have been studied extensively 0- Phase transition has 
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also been observed in a special limiting situation of the 
bus route model that involves two coupled dynamical 
variables one of which, because of the periodic bound- 
ary conditions, is conserved while the other is not. The 
aim of this paper is to investigate the dynamical phases 
and phase transitions in a model of interacting self-driven 
particles 0,0 which has been motivated by the collective 
traffic-like movement of ants on a trail [l(J • 

In this paper we first utilize the close relation be- 
tween the ant-trail model (ATM) and zero-range pro- 
cesses (ZRP) [HII2 to show the existence of a contin- 
uous phase transition from an inhomogeneous jammed 
phase to homogeneous congested phase at a particular 
density in a special limit of this ATM with periodic 
boundary conditions; this phenomenon is argued to be 
closely related to a phase transition in the bus route 
model 0. Then, in order to investigate the boundary- 
induced phase transitions in the ATM we also compute 
its phase diagram with open boundary conditions for a 
wide range of values of the relevant parameters. 

The paper is organized as follows: the ATM is defined 
in Section ^ where the corresponding numerical results 
are summarized briefly. Then we discuss the relation be- 
tween ATM and ZRP in Section ITTT1 where we also give 
analytical calculations based on this analogy, including 
the results in the thermodynamical limit. The phase di- 
agram of the ATM for open boundary conditions is given 
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in Section ITVl Conclusions are drawn in Section Ivl 



II. ANT TRAIL MODEL 

Since our ATM model reduces to the asymmetric sim- 
ple exclusion process (ASEP) in some special limits, we 
begin by defining the ASEP 13]. A fraction of the sites 
of a one-dimensional lattice are occupied initially ran- 
domly by particles each of which can move forward by 
one lattice spacing, with probability q c g, if and only if 
the target site is not already occupied by another par- 
ticle. The updating is done either in par allel or in a 
random-sequential manner. The ASEP [13j with paral- 
lel updating has been used often as an extremely simple 
model of vehicular traffic on single-lane highways Q . 



A. Definition of the ATM model 

Let us now define ATM which is a simple model for an 
unidirectional motion of ants on an existing trail. The 
ants communicate with each other by dropping a chem- 
ical (generically called pheromone) on the substrate as 
they crawl forward [bH Il5l ] . The pheromone sticks to the 
substrate long enough for the other following sniffing ants 
to pick up its smell and follow the trail. The presence of 
pheromone leads to a higher effective speed of the iso- 
lated ants. In ref.8 and 9 we have proposed an extension 
of the ASEP that takes into account this enhancement 
of the effective speed. 

Each site of our one-dimensional ant-trail model rep- 
resents a cell that can accomodate at most one ant at 
a time. The lattice sites are labelled by the index i 
(i = 1,2,...,L); L being the length of the lattice. We 
associate two binary variables Si and o~i with each site 
i where Si takes the value or 1 depending on whether 
the cell is empty or occupied by an ant. Similarly, o~i = 1 
if the cell i contains pheromone; otherwise, Ui = 0. The 
instantaneous state (i.e., the configuration) of the system 
at any time is specified completely by the set ({S}, {cr}). 

Since a unidirectional motion is assumed, ants do not 
move backward. Their forward-hopping probability is 
higher if it smells pheromone ahead of it. The state of 
the system is updated at each time step in two stages. 
In stage I ants are allowed to move. This motion follows 
rules similar to those of the particles in the ASEP 
except that the hopping probability now depends on 
the presence or absence of pheromone at the target 
site. Here the subset {S(t + 1)} at the time step t + 1 
is obtained using the full information ({S(t)}, {<j(t)}) 
at time t. Stage II corresponds to the evaporation of 
pheromone. Here only the subset {o~(t)} is updated 
so that at the end of stage II the new configuration 
{{S(t + 1)}, {a(t + 1)}) at time t + 1 is obtained. In each 
stage the dynamical rules are applied in parallel to all 
ants and pheromones, respectively. 



Stage I: Motion of ants 

An ant in cell i that has an empty cell in front of it, i.e., 
Si(t) — 1 and Si+i(t) = 0, hops forward with 

probability^/^ * = 1' (1) 

(_ q if o- i+1 (t) = 0, w 

where, to be consistent with real ant-trails, we assume 
q<Q. 

Stage II: Evaporation of pheromones 

At each cell i occupied by an ant after stage I a 
pheromone will be created, i.e., 

£7i(i+l) = l if Si{t+l) = l. (2) 

On the other hand, any 'free' pheromone at a site i not 
occupied by an ant will evaporate with the probability / 
per unit time, i.e., if Si(t + 1) = 0, o~i(t) — 1, then 

, , _ J with probability /, , . 

° A + > \ 1 with probability 1 - /. 1 > 

Note that, if periodic boundary conditions are im- 
posed, the dynamics conserves the number N of ants, 
but not the number of pheromones; in that case ATM 
model is a stochastic CA model with two coupled dy- 
namical variables one of which is conserved and the other 
nonconserved. The stationary states of this ATM model 
reduces to those of the ASEP p| with q c g — Q and 
<Zcff = q, respectively, in limiting cases f — and / = 1. 

B. Relation between ATM and the bus route 
model 

The bus route model describes a system of buses 
that move unidirectionally from one bus stop to the next 
on a circular route and, at each bus stop, a bus picks 
up the waiting passengers that arrive stochastically since 
the departure of the last bus from that stop. 

In the bus route model pj the bus stops are represented 
by the sites on a one-dimensional lattice each of which 
may be labeled by an index i (i — 1, 2, L). Two binary 
variables <7j and n are assigned to each cell i: (i) If the 
cell i is occupied by a bus then o~% = 1; otherwise Oi = 0. 
(ii) If cell i has passengers waiting for a bus then r< = 1; 
otherwise r, = 0. Since a cell cannot have simultaneously 
a bus and waiting passengers, a cell cannot have both 
(Tj = 1 and t,- = 1 simultaneously. Each bus is assumed 
to hop from one stop to the next. 

The random- sequential update rules of the model are 
as follows: a cell i is picked up at random. Then, (i) 
if (Tj = and Tj = (i.e, cell i contains neither a bus 
nor waiting passengers), then r% — * 1 with probability A, 
where A is the probability of arrival of passenger(s) at 
the bus stop, (ii) If <Xj = 1 (i.e., there is a bus at the cell 
i) and (Ti-\-i = 0, then the hopping rate \x of the bus is 
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defined as follows: (a) if 73+1 = 0, then fi = ctb but (b) if 
Tj +1 = 1, then p — fa, where a& is the hopping rate of a 
bus onto a stop which has no waiting passengers and (3b 
is the hopping rate onto a stop with waiting passenger (s). 
Generally, (3b < «&, which reflects the fact that a bus has 
to slow down when it has to pick up passengers. One can 
set a.b = 1 without loss of generality. When a bus hops 
onto a stop i with waiting passengers is reset to zero 
as the bus takes all the passengers. Note that, because 
of the periodic boundary conditions, the density of buses 
is a conserved quantity whereas that of the passengers is 
not. 

Note that in the ATM the pheromones "appear" at a 
site when an ant visits the site just as the waiting passen- 
gers "disappear" from a site when a bus visits the site. 
Also note that the "disappearance" of pheromone from a 
site takes place with a probability / per unit time, inde- 
pendent of the ants. Similarly, the passengers "appear" 
at a bus stop with the probability A per unit time in- 
dependent of buses. Therefore, the "appearance" of the 
pheromone is analogous to the "disappearance" of the 
passengers and vice versa. Naturally, we expect the role 
of the parameter / in the ATM to be similar to that of 
A in the bus route model. The presence of passenger is 
analogous to the absence of pheromone. Thus, the ATM 
is equivalent to the BRM with parallel updating 

Interestingly, the queueing of the buses in the bus 
route model, in turn, is very similar to the bunching 
of macrosteps during the growth of faceted crystals be- 
cause of impurity absorption |17|. If y n (t) denotes the 
position of the particle (representing a bus or an ant or a 
macrostep) at time t, and r n is the time elapsed since the 
particle ahead (which is labelled by n+1) passed the same 
position y n (t), one can write [13 y n (t) = y n +i(t - T n ) 
and corresponding speed v n (t) = V{T n (t)) is a function 

of T„(t). 



C. Numerical results 

We shall use the symbols L and M to denote the sys- 
tem size and the number of ants, respectively; hence 
p = M/L is the density of ants on the trail. The most 
important quantity of interest in the context of flow prop- 
erties of the traffic models is the fundamental diagram, 
i.e., the flux-versus-density relation, where flux F is the 
product of the density p and the average speed v. Be- 
cause of the similarity of the ATM with traffic models, it 
would be interesting to draw the fundamental diagram of 
ATM and compare the traffic-like collective movements 
of ants with vehicular traffic. 

The fundamental diagrams of ATM is given in our ear- 
lier papers (see, for example, Fig. 2(b) of ref.8 or Fig. 2(b) 
of ref.9). First of all, the diagram does not possess 
particle-hole symmetry for any / in the range < / < 1; 
the particle-hole symmetry observed in the ATM in the 
special cases of / = and / = 1 is a consequence of the 
fact that in the two special cases / = and / = 1, as 



pointed out earlier, the ATM becomes identical to the 
ASEP, with parallel updating, corresponding to the ef- 
fective hopping probabilities q c s = Q and g e ff = Q, re- 
spectively. 

Most important feature of the fundamental diagram 
is that, over a range of small values of / (/ < 0.01), it 
exhibits an anomalous behaviour; this is a consequence 
of the fact that, unlike common vehicular traffic, the av- 
erage speed of the ants in the ATM is a non-monotonic 
function of the density of ants on the trail (see, for ex- 
ample, Fig.3(b) of ref.8 or Fig.2(a) of ref.9). 

The ATM model also exhibits interesting coarsening of 
clusters of ants starting from random initial conditions 
(Fig.l). 



x 




FIG. 1: Spatial-temporal behaviours of loose clusters in the 
low density case (p = 0.16). Parameters are Q — 0.75, q = 
0.25,/ = 0.005). We see the loose clusters emerge form the 
random initial configuration, and will eventually merge into 
one big loose cluster if we take sufficiently long time. 



III. PHASE TRANSITION IN ATM WITH 
PERIODIC BOUNDARY CONDITIONS 

In ref.9 we have developed a formalism by introduc- 
ing the loose cluster approximation and successfully cap- 
tured the non-monotonicity heuristically with the help 
of analytical results of ASEP, which is an exactly solv- 
able stochastic model. In this paper, we study the phe- 
nomenon in detail by utilizing the analogy of ATM with 
the zero range process (ZRP), which is exactly solvable 
[TTL IT^ . The ZRP is a particle-hopping model where 
the hopping probabilities do not depend on the state of 
occupation of the target site. 

A. Analytical results for finite systems: ZRP and 
ATM 

We now first show the relation between the ATM and 
the ZRP [HII3- Note that > in thc A ™ the effective 
hopping probability u of the can be exactly expressed as 

Si 

u = q(l - g) + Qg, (4) 
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where g is the probability that there is a surviving 
pheromonc on the first empty site of a gap in front of 
an ant. Obviously, g(t) depends on the time t elapsed 
since the ant immediately in front left the target site. It 
is straightforward to see that 



S(* + 1) = (1 -/)<?(*) 



(5) 



Alternatively we can express the probability g(t) as a 
function g{x) of the gap x between the ant under consid- 
eration and that immediately in front of it. Recall that 
the average speed of the ants is v. Since the time inter- 
val between the passage of successive ants through any 
arbitrary site is x/v, and since equation JJjJ holds for all 
t, we obtain 



g(x) = (1 - /)*/« 



(6) 



after iterating the mapping J5J x/v times. 

Thus, the ATM can be regarded as an ASEP involv- 
ing fictitious particles where, u(x), the effective hopping 
probability of a particle having a gap x in front, depends 
on x through the relation Q 



u(x) = q+(Q-q)(l-fr/\ 



(7) 



Note that in the case of the original ASEP, u(x) — con- 
stant for x > and u(0) = 0. The form Q) is approxi- 
mate description of the hoppings of the ants in the ATM, 
but all the results we derive for this form of u(x) are ex- 
act. 

As long as / remains non-zero, u(x) decreases with 
increasing a; at a rate much faster than that required 
for a phase transition to occur 18]. However, since v 
depends on / and x, the trend of variation of u{x) with 
x is difficult to infer without going through a detailed 
analysis which we present below. 

The configurations of this ASEP can be uniquely de- 
scribed by the gap configurations {x\,X2, ...,xm}- Note 
that in this ASEP the phcromones do not appear explic- 
itly, but their effect enters through the /-dependence of 
u(x). It is well known ^1 that any ASEP, where the 
hopping probability of the particles u(x) depends on the 
gap x in front of the particle, can be mapped onto the 
ZRP. The advantage of mapping the ATM onto the cor- 
responding ZRP (with parallel updating) is that the the 
stationary state of the ZRP is given by product measure 
in spite of the inter-particle interactions. 

In order to avoid any possibility of confusion, we now 
point out the commonly used notation for ZRP |l9j |. In 
that notation for ZRP, the lattice consists of M sites each 
of which can hold an integer number of indistinguishable 
particles. The configuration of the system is specified 
by the occupation numbers {ni, n,2, ■-, Um}- The total 
number of particles is L. On mapping this ZRP onto 
the ASEP, each site of ZRP becomes a particle in ASEP 
while the total number of sites of the ASEP becomes 
identical to the total number of particles in ZRP. There- 
fore, the expressions of some quantities like, for example, 



the average speed may appear slightly different in the 
two notations, but are completely equivalent. 

Following the treatment of the bus route model in ref.7, 
now the steady-state probability PmUx^}) of finding the 
ATM in a gap configuration {x\,X2, ...,%} is given by 
a product of factors h(x^), i.e., 



p (s w n^Lx %m) 



(8) 



where the partition function Z(L, M) is given by the re- 
lation 



M 



Z(L,M)= H h M5[E^ -(L-M) 

(9) 

The Kronecker detla in @ has the standard meaning, 
i.e., S(a,b) = 1 if a = b and 6(a,b) — otherwise. Note 
that Z(L,M) is just a normalization which ensures that 
the sum of the probabilities Pm({^}) over all possible 
gap configurations {x\,X2, ...,Xm} is unity. 

The single-gap probability p{x), i.e., the probability 
that there is a gap of size x in front of an ant, is obtained 
from the M-gap probability P({Xfj,}) by using 



p(x) = Pi (x) 



Y S[x 2 + ... + x M -{L- M - x) ) h{x) Yl Kx^) 

fi=2 



h(x 



Z{L, M) 
Z(L-x - 1,M- 1) 
ZjL, M) 



(10) 



where Z(L — x — 1, M — 1) is the partition function for a 
system from which the first site has been removed. 

Note that Z(x, 1) = h(x-l) and Z(x, x) = h(0). Then, 
in principle, if h{x) is known, the partition function Z can 
be obtained by using the recursion relation 

L—M 

Z[L, M) = ^ Z(L — x — 1, M — l)h(x), (11) 

The ATM is formulated with parallel update. The cor- 
responding form of h(x) is known to be given by |20| 



l-u(l) 

h{x) = { !-"(!) -pr 1 - u(y) 

1 ~ M ( a; ) yJi u ^ 



for x = 0, 
for x > 



(12) 

Our aim is to calculate the average velocity v of ants 
given by 



L — M 



u(x)p(x). 



(13) 



In order to compute v we need to compute p(x) and u(x). 
However, the computation of p(x) requires h(x) which, in 
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turn, depends on u(x). On the other hand, u{x) is given 
by equation J7J which involves average speed v. There- 
fore, in order to compute v self-consistently, we begin 
with the initial approximate estimate v — q and, hence, 
compute u(x) and p(x) to get better estimate of aver- 
age speed v from equation Through this iterative 
process, we estimate v self-consistently and, hence, the 
fundamental diagram. 

Fundamental diagrams are given in Fig. 2 with L = 100 
and L = 200. The data points, shown with various 
symbols in Fig. 2 and connected by dashed curves, have 
been obtained from computer simulations of the ATM, 
whereas the continuous curves are the corresponding the- 
oretical predictions. We find that the theoretical curves 
are almost identical to the numerical ones, thus confirm- 
ing that the ZRP successfully describes the steady state 
of the ATM. 
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FIG. 2: Comparison of the theoretically obtained funda- 
mental diagrams of the ATM with for the system sizes (a) 
L = 100 (continuous curve in (a)) and (b) L = 200 (con- 
tinuous curve in (b)) with the corresponding numerical data 
obtained from computer simulations (represented by the sym- 
bol +). The numerical data obtained from simulations of 
systems of size L — 1000 are also plotted (represented by 
the symbol x) in both (a) and (b) to show the trend of 
variation with increasing L. The common parameters are 
Q = 0.75,g = 0.25,/ = 0.005. 



B. Thermodynamic limit of ATM 

Next we discuss the thermodynamic limit of the ATM, 
that is, the case L — > oo. From Fig. 2 we see that 
the curve shows sharp increase near the density region 
0.4 < p < 0.5, and the tendency is expected to increase 
with the increase of L. Thus it is important to study 
whether there is a (second order) phase transition from 
an inhomogeneous jammed phase, in which the ants are 
bunched together, to homogeneous congested phase with 
the increase of the density of the ants in the thermody- 
namic limit. 

Using the integral representation 



5(ji,v) 



ds 



2iri s v+1 



(14) 



of the Kronecker delta, we rewrite the partition function 



Z( L, M ) as 



Z{L,M) 



ds 



G{s) 



M 



2-nis \sVp-i 

where G(s) is the generating function of h defined by 



(15) 



G(s) = J2h(x)s*. 



(16) 



We evaluate the integral in (|15|) in the L — > oo limit, keep- 
ing M/L = p as constant, by the saddle point method. 
The partition function reduces to the form 

Z(L,M) ~ exp(MmG(z) - {L - M)\nz) (17) 

where the saddle point s = z is given by 



1 

P 



1 = z 



d\nG{z) 

dz 



(18) 



The equation l|18|) may also be regarded as the relation 
that defines z. 

Using the approximate form (|17fl of the partition func- 
tion, we have 



P ( X ) ~ m z{L ~ x ~ 1,M ~ l) ~ ^ 

Then, substituting (|19|l in expression (|13|l we get 

u(x)h(x) 



G(z) 



(19) 



(20) 



for the average speed v in the thermodynamic limit. 

Next we study the properties of the generating func- 
tion G(z) in detail to examine the possibility of phase 
transition of the ATM. Since lim^^oo u(x) — q for all 
/ > 0, the condition that G converges is given by 



h(x + l)z 1 
urn sup — — — — < 1 . 

x^oo h{x) 



(21) 



Using i|12|) for h(x) the convergence condition reduces to 



h(x + l)z 1 - q 
lrni sup — — — — = z < 1. 

x^oc h[x) q 

Thus, G converges in the range 

q 



< z < z c 



1-9' 



(22) 



(23) 



Using and (|23|) in l|l(j[l we find that the critical value 
of G at z = z c is given by 



1 — u(l) -rj / 1 — u(y) q 



x=l 



0) yJi V u (v) 1 - ? 



(24) 
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For all / > 0, this sum diverges for the class of u(x) 
which decays to q more rapidly than the function (1 + 
c/x)q, where c > 1 — q. In the ATM, u(x) — q + (Q — 
q)(l — f) x / v , which decays exponentially to q as x — >• oo. 
Hence we conclude that there is no phase transition in 
the ATM for / > 0. This is because, from lfT%)l. we 
have p = 1 when z = 0, and p = at z = z c if G{z c ) 
diverges. Thus, as long as f > 0, in the entire density 
region < p < 1 there is no singularity in G and, hence, 
no phase transition in the ATM. 

The situation is drastically different in the limit / — » 0. 
In this limit, u(x) = Q, and then G{z c ) approaches the 
finite limit 



Urn G{z) = 

z — *Zc 

Thus, in the limit / - 
at the critical density 



Q(l-q) 



Q. (25) 



Q(l - q) - q{l - Q) 
-> a phase transition takes place 



Q-q 



(26) 



which is obtained from l|18J). 

In order to get insight into the nature of the phases 
involved, we now calculate the average speed v = v c at 
z = z c . Since u(x) — * Q as / — ■> 0, h(x) — > (1 — Q)/Q 
and, hence, from l{2Tjjl we get 



E 



tG(z c )V 



(27) 



Substituting (|25|l into l|27|) . we obtain v c = q. 

At first sight the two facts, namely, u{x) = Q and 
v c = q may appear mutually contradictory. But, the 
consistency of these is a consequence of the fact that the 
mutual hindrance leads to the smaller average speed v c — 
q in spite of the higher hopping probability u{x) = Q. 

It should be noted j§| that lj2rj|) is also obtained by the 
intersection point of the line F — v c p and the curve 



F=-(l-y/l-4Qp(l-p)) 



(28) 



that describes the fundamental diagram of the ASEP 
with parallel updating. Also note that the limits L — > oo 
and / — > do not commuteQ- If we take / — > be- 
fore L — > oo, then we apparently have l|28|l . which cor- 
responds to the situations in our numerical simulations. 
On the other hand, if we take / — > after L — ► oo, then 
we have the thick curve in Fig. 3. Thus, in the latter 
case, we have proved that the anomalous variation of the 
average velocity with the density disappears. 
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FIG. 3: The thick curve (with no decorating symbols) repre- 
sents the fundamental diagram of the ATM corresponding to 
/ — » after taking the thermodynamic limit L — * oo. Note 
that these two limits do not commute. The simulation data 
corresponding to / = 0.0005(0), 0.001(o), 0.005(»), 0.01(A), 
0.05(a), 0.10(x), 0.25(+), 0.50(*). 



boundary conditions are more realistic. Therefore, in 
this section we study the phases and phase transitions 
in the ATM with open boundary conditions. 

Suppose a and (3 denote the probabilities of incoming 
and outgoing particles at the open boundaries per unit 
time. The phase diagram of the ASEP in the a — /3-planc 
has been investigated exhaustively 2] . In this section we 
report the effects of varying the pheromone evaporation 
probability / on the phase diagram of the ASEP with 
parallel updating and open boundary conditions. 
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FIG. 4: The phase diagram of the ATM with open boundary 
conditions in the a — /3-plane for several different values of 
the pheromone evaporation probability / (0 < / < 1). The 
values of the hopping parameters are Q = 0.75, q = 0.25. 



IV. PHASES OF ATM WITH OPEN 
BOUNDARY CONDITIONS 

So far we have considered the ATM with only periodic 
boundary conditions. However, for ant-trails the open 



Just as in the case of the ASEP, we found three dif- 
ferent phases, namely, the high-density phase, the low- 
density phase and the maximal current phase of the ATM 
for all / (see Fig. 4). The line separating the low-density 
phase and the maximal current phase is given by a c (f) 
whereas that corresponding to the separation between 
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the high-density phase and the maximal current phase is 
given by (3 c (f). 
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FIG. 5: The critical rates, a c , (denoted by +), which have 
been obtained from computer simulations of the ATM with 
open boundary conditions, are plotted against the pheromone 
evaporation probability /. The discrete data points denoted 
by the symbols x and * denote the average speeds t>i/2(/) 
obtained from computer simulations of the same model but 
with periodic boundary conditions for system sizes L = 100 
and L = 500, respectively. The corresponding theoretical 
predictions for «i/2(/) of our ZRP-based theory of ATM are 
plotted for the system sizes L — 100 (the continuous curve) 
and L = 500 (the dashed curve), respectively. The values of 
the hopping parameters are Q = 0.75, q — 0.25. 



The variation of the critical point a c (f) = (3 c (f) with / 
is shown in Fig. 4. In order to understand this observation 
let us first examine the limiting balues at / = and 
/ = 1. It is well known |2lJ that for ASEP with hopping 
probability q e e, 



= e = i- v/T - 



q C ff- 



Therefore, we must have 



a c (f = 0) = c (f = 0) = 1 - y/l-Q 
a c {f = \)=p c {f = l) = l-y/T=q~ ; 



(29) 



(30) 



these are consistent with the numerical data a c (f = 0) = 
0.5 and a c (f = 1) = 0.133 for Q = 0.75, q = 0.25 shown 
in Fig. 4. 

Next we try to understand the detailed variation of 
&c = (3 C with / by utilizing the relation between the ATM 
model with periodic and open boundary conditions. Note 
that, in the maximal current phase in the ASEP with 
open boundary conditions and parallel updating |2l| the 
current is given by 



J = 



1 - V 71 - <hs 



(31) 



(in ASEP with parallel updating) while the correspond- 
ing bulk density is given by |2l| 



and, consequently, 
should be 



the corresponding average speed 



l-vr 



9cff = a c 



(33) 



In order to check the validity of this argument, we 
have computed the average speed V1/2 corresponding to 
p = 1/2 by (a) computer simulation of the ATM with pe- 
riodic boundary conditions and (b) using the ZRP-based 
theory, mentioned in the preceeding section. All these 
data are plotted in Fig. 5. The simulation data for u x / 2 (/) 
are in good agreement with the simulation data for a c . 
But, there are significant differences between these data 
and Vi/2(f) obtained from our ZRP-based theory in the 
small / regime. We believe that this discrepancy arises 
from the boundary effect. If periodic boundary condi- 
tions are imposed, at low densities, the leading ant in the 
loose cluster can smell the pheromone which is dropped 
by the last ant in the same cluster. However, this effect 
would disappear when the periodic boundary conditions 
are replaced by open boundary conditions. 

Since we have numerically estimated a c in the ATM 
as a function of / by carrying out computer simulations, 
we utilized the relation (|2*5|l to get the /-dependence of 
the effective hopping probability q c s(f) in the ATM (see 
Fig. 6). The two limiting values lim/_>n Qes(f) = Q and 
lim/_>i <7cff (/) = 9 as well as the nature of the variation 
of q e g with / are fully consistent with one's expectation 
based of physical arguments mentioned earlier in this pa- 
per. 
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p(x = L/2) = 1/2 



(32) 
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FIG. 6: /-dependence of the effective hopping probability, 
extracted by using the equation 1291 from the numerical 
data obtained by computer simulation of the ATM with open 
boundary conditions. 



V. CONCLUDING DISCUSSIONS 

The close similarities between phase transitions in non- 
living systems and that of foraging behaviour on the ant- 
trail was pointed out by Beckman et al.JU- However, 
their study was concerned with the nature of foraging, 
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namely, the possibility of a transition from a disordered 
foraging behaviour (i.e., foraging without a pheromone 
trail) to ordered foraging (i.e., trail-based foraging). In 
contrast to this problem of the formation of the trail pat- 
tern |23| . we have studied the phases and phase transi- 
tions exhibited by the system when the collective traffic- 
like movement of ants take place on an existing trail. 

The dynamical phases and non-equilibrium phase tran- 
sitions in systems of interacting self-propelled parti- 
cles have been among the most challenging problems 
of investigation in non-equilibrium statistical mechan- 
ics. In contrast to equilibrium systems, these intrinsi- 
cally non-equilibrium systems exhibit phase transitions 
even in one-dimensions with only short range interac- 
tions. In this paper we have studied the phases of a one- 
dimensional model motivated by the collective traffic-like 
movements of ants on an existing trail 0, 0- In this 
model the ants are represented by self-propelled particles 
which, in addition to the hard-core repulsion, interact 
indirectly via pheromone. The model is, thus, charac- 
terized by two coupled dynamical variables, representing 
the ants and the pheromone. 

In our earlier works |8|,|9j we had shown that the homo- 
geneous mean-field approximations cannot capture the 
non-monotonic variation of the average speed with the 
density of the ants. Even the loose-cluster approxima- 
tion that we developed|9j could account for the simula- 
tion data with limited accuracy. In this paper we have 
reported our new quantitative results on this model with 
periodic boundary conditions; these results have been de- 
rived by utilizing the analogy with ZRP. Moreover we 



have shown that there is a phase transition in the ther- 
modynamic limit in this model, albeit in a special limit 

In our earlier published works 0, Q we imposed peri- 
odic boundary conditions. However, in order to capture 
the real ant-trails in nature imposition of open boundary 
conditions seems more appropriate. In equilibrium sta- 
tistical mechanics, the boundary conditions do not play 
any role in the phase transitions which, strictly speak- 
ing, take place only in the thermodynamic limit. In con- 
trast, boundary- induced phase transition [(| is a well- 
known phenomenon in non-equilibrium statistical me- 
chanics. The phase diagrams of the ASEP with open 
boundary conditions is now well established. Since the 
ATM can be regarded as an extension of the ASEP, our 
investigation of the phase diagram of the ATM reported 
here illuminates the effects of varying the pheromone 
evaporation rate / on the phase diagram of the ASEP. 
In this paper we have focussed attention on phase transi- 
tions in an idealized single-lane ATM. The phase transi- 
tions in a more elaborate two-lane ATM will be reported 
later j^. It would be interesting to test the predictions 
of the model, particularly the non-monotonic variation 
of average speed with density, by repeating experiments 
of the type reported in ref.10 to get sufficiently accurate 
data for real ants. 
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